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Abstract 



The calculation of the lifetime of quasiparticles in a QED plasma at high tem- 
perature remains plagued with infrared divergences, even after one has taken into 
account the screening corrections. The physical processes responsible for these di- 
vergences are the collisions involving the exchange of very soft, unscreened, mag- 
netic photons, whose contribution is enhanced by the thermal Bose-Eisntein oc- 
cupation factor. The self energy diagrams which diverge in perturbation theory 
contain no internal fermion loops, but an arbitrary number of internal magneto- 
static photon lines. By generalizing the Bloch-Nordsieck model at finite tempera- 
ture, we can resum all the singular contributions of such diagrams, and obtain the 
correct long time behaviour of the retarded fermion propagator in the hot QED 
plasma: Sji(t) ~ exp{— aTt lniv p t}, where lo p = eT/3 is the plasma frequency and 
a = e 2 /47r. 
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It is commonly assumed that the notion of quasiparticles is a useful concept in 
the description of ultrarelativistic plasmas. This implies that the damping rate 7 of 
quasiparticle excitations, obtained from the exponential decay of the retarded propaga- 
tor, Sn(t, p) ~ e -Mp)* e -7(p)* ? j s sma ll compared to the quasiparticle energy w(p). In 
hot gauge theories, the typical energy of the relevant quasiparticles is the temperature 
T, while one expects 7 ~ g 2 T [l], [|, ||, where g is the gauge coupling (in QED, g = e 
is the electric charge). The same damping rate is expected for the collective excitations, 
whose typical energies are ~ gT p. This suggests that, indeed, in the weak coupling 
regime, quasiparticles are well defined, and collective modes are weakly damped. How- 
ever, the computation of 7 in perturbation theory is plagued with infrared divergences, 
which casts doubt on the validity of these statements [1-7]. We shall show in this pa- 
per that the infrared divergences which remain after the screening corrections have been 
taken into account are due to collisions involving the exchange of very soft, unscreened, 
magnetic photons. Such divergences occur in all orders of perturbation theory and can 
be eliminated only by a non perturbative calculation of the quasiparticle propagator. 
We present in this letter such a non perturbative analysis, based on an extension of the 
Bloch-Nordsieck model at finite temperature. The final result is that quasiparticles do 
indeed exist, although they do not correspond to the usual exponential decay indicated 
above, but to a more complicated behaviour, £r(£, p) ~ e -iw(p)*e -aT * lnw p* ) w h ere u ^ j s 
the plasma frequency. 

We shall consider in most of this letter the case of hard fermions, with momenta 
p ~ T, and restrict ourselves to QED. (The case of collective excitations will be dis- 
cussed briefly at the end.) In perturbation theory, the damping rate is obtained from the 
imaginary part of the self-energy on the unperturbed mass-shell uo = p : 



To bare one-loop order, 7 = by kinematics : the on-shell fermion can emit (or absorb) 
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only a space-like photon (q 2 = q 2 — q 2 < 0), and there are no such photons at tree level 
(a more complete discussion of this point will be presented in a longer publication ||). 
However, through the Landau damping mechanism, the spectral function of a photon in 
the plasma does acquire a non vanishing weight for space-like momenta f9j. One takes 
this into account by resumming the so-called hard thermal loops (HTL) in the internal 
photon line of the one loop self-energy, and obtains in this way the leading contribution to 
the damping rate [l], |2|. Physically, this contribution represents the effect of the collisions 



induced by one photon exchange [10]. Note that the corresponding collision rate would 



be quadratically divergent without the screening corrections [1G, 11]. Such corrections are 



sufficient to render finite the transport coefficients [10], or the collisional energy loss [3 I 



but not the total interaction rate []12 |, nor the damping rate [1-7]. These quantities, except 
for the damping rate of a quasiparticle with zero momentum 0, [13|] , remain logarithmically 
divergent, both in abelian and non-abelian gauge theories. In QCD, this problem is 
commonly bypassed by advocating the IR cut-off provided by a possible magnetic mass 
~ g 2 T (see e.g. H). But such a solution cannot apply for QED where one expects no 



magnetic screening ||14|| . 



In the imaginary time formalism, and in the Coulomb gauge (the one-loop result 
for 7 is gauge independent |j| ) , the resummed one loop diagram is evaluated as H 



= -9 2 T £ J^- 3l ,S (p + q) lu *D^(q) 



(2nY 

2 f d 3 q [+°°dk [+™dq „ N (go) + w(fc ) 

9 / 7^3 / / 77-Po(^)7 M f7^ P M (<?) -r • (2) 

J {Ztt) 6 J-oo 2tt J-oo In k — q — p 



In this equation, k = p + q, p = iuj n = i(2n + l)vrT, and uj m = 2TrmT, with integers 
n and m. We denote by po(k Q ,k) the spectral density of the (hard) fermion, and by 
*Pi,t(qo,q) those of the (soft) longitudinal (I) and transverse (t) photons. According to 
eq. (|IJ), one obtains the damping rate 7 from the imaginary part of *£, after analytical 
continuation of po to the real axis (po —* w + iv)- O ne thus gets 7(0? ~ p) = g 2 p^(u,p), 
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where (cos 9 = p • q) 



~N(q )+n(k )]p (k) + (1 - cos 2 0)*pt(g)) , (3) 

is the thermal phase-space for the decay of the fermion. The infrared singular region is 
that of small photon energy-momentum go, q <C gT. Since u, p ~ T, in studying this 
region, one can expand |p + q| ~ p + qcos9 and N(q ) + n(k ) ~ N(q ) ~ T/q ^> 1. 
One obtains then 



where the 5-function has been used to perform the angular integration, and 8E = \uj — 
p\ > 0. Note the restriction to space-like momenta, |go| < 9 : only the off-shell photons 
contribute to the fermion damping, which is consistent with the physical interpretation of 
7 in terms of collisional processes, as alluded to earlier. The evaluation of the go-integral 
can be done easily by using sum rules satisfied by the spectral function ]5). The singular 
contribution comes from the transverse component, for which we have 

where u> t (q) and z t (q) are respectively the position and the residue of the transverse 
plasmon pole. As q — > 0, z t (q)/uf(q) — > l/cu 2 , so that for q <C u p ~ gT the integral 
is dominated by the term 1/g 2 . In fact, the function *Pt(qo,q)/qo is strongly peaked at 
g = 0, and in the calculation of the integral (^) it can be replaced by the following 
approximate expression: 

— Pt{q < q ) = — 6 ; 7z 2 ,.. 2 — > q -*o — 5{q ). 6 

g 2 q b + (37ru^g /4) 2 q l 

Keeping only this leading term, and ignoring the non-singular contribution of the electric 



modes, one finds 

1 {5E)=g 2 p<$>{5E)^ y — dqq6(q-5E)- = aT\n-^, (7) 

47r jo q oti 

with a = g 2 /4:TT. Thus, the singularity of 7 when SE = 0, i.e., when the fermion is 

on-shell, can be attributed to the low frequency magnetic photons produced by Landau 

damping, whose contribution is enhanced by the Bose-Einstein occupation factor T/q . 

Note that since the energy conservation requires qo = q cos 9, these photons are emitted, 

or absorbed, at nearly 90 degrees. 

The upper cut-off u p ~ gT in eq. (^) accounts approximately for the terms which 
have been neglected when keeping only the 1/q 2 contribution to the r.h.s of eq. (|5]): as 
q ^> Up, the r.h.s. of eq. @ is decreasing like Wp/g 4 , so that the exact integrand in eq. @ 
is indeed cut-off at q ~ u p . Note that, as long as we are interested only in the coefficient 
of the logarithm, the precise value of this cut-off is unimportant. The scale uo p however 
is uniquely determined by the physical process responsible for the existence of space like 
photons, i.e., the Landau damping. This is important since, as we shall see later, it is the 
scale which fixes the long time behavior of the retarded propagator. 

In the forthcoming discussion, we shall take advantage of a major technical simplifi- 
cation which is based on the following observation: the logarithm in eq. (|7|) arises entirely 
from the magnetic contribution of the term q = iuo m = in the Matsubara sum of eq. (Q) 
f|], . The analytic continuation of this term to real energy ( po —* 00 + if]) is well defined, 
and it yields the following imaginary part 

1 r d 3 o 1 1 

--trOdm-EoW) . -V^/^ 5? Im ( ^ + . i;)2 _ (p + q)2 

47T Jo q 

where the approximate equality means that only regular terms have been dropped. As 
anticipated, this coincides with eq. (0). Note that, since the static transverse propagator 
is not modified by the HTL corrections, eq. (|j) could have been obtained directly from 

4 



the term q = iuo m = in the bare one-loop self-energy, provided one introduces the same 
upper cut-off u p in the (/-integral. 

Some remarks are needed at this stage. At the bare one- loop level, 7 = 0, as 
mentioned. Thus, in the bare one loop calculation, there is a cancellation between the 
logarithmic divergence of the static mode, and that of the sum of all the non-static ones 
||. The inclusion of the screening corrections on the photon line removes the divergence 
of the contribution of the nonstatic modes, leaving a divergence in the static contribution 
only. When ui — > p, the integrand in eq. (H) becomes proportional to 5 (cos 6), so that 
the singular contribution of the static mode originates from the emission or absorption 
of static magnetic photons at 90 degrees. As we have seen, these processes represent the 
dominant contribution of the resummed one loop calculation. That they can be taken 
into account by what looks like a bare one loop calculation is due to the fact that the 
integrated transverse photon spectral weight (||) is 1/q 2 at small q, i.e., is indistinguishable 
from the propagator of a static photon. Finally, as already mentioned, the source of the 
space-like photons is Landau damping; this is the physical origin of the cut-off u p to be 
imposed on the one loop integral of eq. 

Power counting indicates that mass-shell divergences do occur in higher orders as 
well. The most singular diagrams are those where the fermion propagator is dressed 
by static magnetic photon lines. When evaluated on the tree-level mass-shell u = p, 
and in the presence of an IR regulator //, such diagrams generate contributions of the 
type (g T/fi) n , where n is the number of loops. Such power-like IR divergences are 
analogous to those identified in the analysis of the corrections to the screening mass in 



T6j, and their presence signals a breakdown of perturbation theory. The calculation of 
the fermion propagator in the IR singular domain requires therefore a non perturbative 
approach, to which we now turn. 

Since internal fermion loops are not essential anymore at this stage, we may restrict 
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ourselves to the "quenched approximation" involving only magnetostatic photons EBI. In 
this approximation, the imaginary-time fermion propagator can be written as the following 
functional integral 

S(x,y) = Z - 1 |[dA]G(x,y|A)exp|-i (A, D^A)^ , (9) 

where G(x,y\A) is the fermion propagator in the presence of a background gauge field, 
and 

(a,d 1 a) o = iy dW^(x)Z^(x-y)^(y). (10) 

In this equation, the factor 1/T has its origin in the restriction to the zero Matsubara 
frequency. The propagator D 0i j is that of a free static photon, 

Di j (x) = J (05 eiq ' X ^'(q) (11) 

with 

Dl j (q)= ^ + (A-l)£f, (12) 

gZ gl 

in an arbitrary Coulomb or covariant gauge (A = corresponds to both the Landau and 
the strict Coulomb gauges). 

The fermion propagator in a static background field depends only on the time dif- 
ference Xq — i/q, G(x,y\A) = G(xq — y , x, y| A). Its Fourier transform can be analytically 
continued in the complex energy plane, and the resulting function coincides, in the up- 
per half plane, with the retarded propagator. It is then convenient to take the Fourier 
transform of eq. ffl). To this aim, we write (p^ = (itu n , p), u n — (2n + \)-kT): 



G(x,y\A) = TY: J ^e-*-^G?(p;x|A), (13) 
and similarly for S. We obtain 

S(p) = Z, 1 J [dA] G(p; x| A) exp |-i (A, D^A)^ . (14) 



Since the energy po enters eq. (|T4| ) as an external parameter, the continuation to real ex- 
ternal energy po — > u + irj, and the Fourier transform to real time, can both be performed 
before doing the functional integration |S|]. This offers the possibility to calculate the re- 
tarded propagator Sn(t, x) by directly inserting an approximate expression for Gr(x, y\A) 
in the functional integral (H). 

In the kinematical regime of interest, an approximate expression for Gr{x, y\A) is 
obtained by neglecting the recoil of the fermion in the successive emissions or absorptions 
of very soft photons. More precisely, we note that inside the diagrams generated by 
eq. (H), we can approximate the fermion propagators by 

c / . x -^7o + (p + q) • 7 - 1 70-P-7 , 1K v 
5 (w,p + q) = — ( — 2 - —■ n ' 15 

where q is a linear combination of the internal photons momenta and v = de p /dp (v = 
p for the ultrarelativistic fermion). This is the familiar structure encountered in most 
treatments of IR divergences in QED (see, e.g., [0]) and which is economically exploited 
within the Bloch-Nordsieck model (see, e.g. |I8|). In this model, the propagator with the 



desired IR structure is obtained in coordinate space as the solution to 

- i(v ■ d x )G (x,y) = 8(x-y), (16) 

where = (l,v). In the presence of an arbitrary background field, the corresponding 
propagator, G(x,y\A), satisfies the following equation (D^ = + igA^) 

-i(v-D x )G(x,y\A) = 5(x-y), (17) 

which can be solved exactly. For retarded boundary conditions, and for static fields: 

G R (x,y\A) = t 6(x -y )5W(x-y-v(x -y ))U(x,y) (18) 

where U(x, y) is the parallel transporter 

U(x, x — vt) — exp jig J ds v • A(x — vs)j . (19) 



The retarded propagator Sr(x — y) is calculated by inserting the expressions (|IBD, 
(19) of Gr(x, y\A) in the functional integral (|i~4"|). It can be written as 

S R (t,p) = t6(t)e~ u ^-^A(t), (20) 

where 

A(t) = Z, 1 J [dA] U(x, x - vt) exp j-i (A, D 1 A)J (21) 

contains all the non-trivial time dependence. The functional integral is easily done: 

A(t)=exp|^T jf t ds 1 jr*ds 2 ^DS''(v( S2 -3 2 ))^|. (22) 

In this equation, Dq (x) is the coordinate space representation of the magnetostatic photon 
propagator (see eq ([TT|)). 

The Si and s 2 integrations in eq. (|22|) can be done by going to the Fourier represen- 
tation: 

A(f) = exp |-, 2 T / ^3^5 (l " cost(v ■ q))} , (23) 

where -Dg^q)?;- 7 = -D(q). The integral in eq. (|23| ) is identical to that one would get in 
the Bloch-Nordsieck model in 3 dimensions. It has no infrared divergence, but one can 
verify that the expansion of A(t) in powers of g 2 generates the most singular pieces of the 
usual perturbative expansion for the self-energy ||. One can also verify that the integral 
in eq. ( p3"D presents an ultraviolet logarithmic divergence. However, one should recall that 
the restriction to the static photon mode implies that such an integral is to be cut off at 
momenta q ~ u p (cf. the discussion after eq. (0)). 

The calculation of A(t) is most simply done using the coordinate space representa- 
tion (^2|). In the Feynman gauge A = 1, we use v l Dq (x) = l/Airx and obtain 



~ f t ds 1 [ t ds 2 v i D i j (v(s 2 -s 2 ))v j 
2 Jo Jo 



1 f l , r\ #Oi - s 2\ - i/w p ) t , 

— / d Sl / ds 2 11 21 1 p> ~ — \nu p t. 24 
37T JO JO Si — s 2 Air 



where the ultraviolet cut-off is introduced in the function 9(\si — s 2 | — 1/uj p ). This ex- 
pression is valid for times t ^> l/u p . It is insensitive to the specific procedure which is 
used to implement the cut-off || , and independent of gauge fixing (the gauge-dependent 
contribution to the integral above, as given by the last term of the photon propagator 
(0), is (A — 1) g^r, and therefore sub leading). 

At times t ^> l/u p the function A(t) is thus of the form 

A(uj p t > 1) ~ exp(-aTtlntUpt). (25) 

A measure of the decay time r is given by 

- = aT\nu p r = aT (ln^ - Inln^ + ...V (26) 
r \ aT aT ) 

Since aT ~ guj p , r ~ 1 / (g 2 T ln(l / g)) . This corresponds to a damping rate 7 ~ 1/r ~ 
g 2 Tln(l/g), similar to that obtained in a one loop calculation with an IR cut-off ~ g 2 T 
(cf. eq. ©). 

However, contrary to what perturbation theory predicts, A(t) is decreasing faster 
than any exponential. It follows that the Fourier transform 

/oo roo 
dte- wt S R {t,p) = i / dte tt{uJ - v - p+iri) A(t), (27) 
-00 Jo 

exists for any complex (and finite) uj. Thus, the retarded propagator Sr{oj) is an entire 
function, with sole singularity at Imw — > —00. 

The previous analysis can be extended to quasiparticles with soft momenta, 
p ~ gT, for which the same infrared difficulty arises H]. Because of the collective nature 
of the soft quasiparticles, the HTL resummations to be performed in this case are more 
elaborate 0. It can be shown however that the leading divergences are again confined 
to the magnetostatic photon sector 0, and occur in diagrams which involve only the 
3-point photon-fermion vertex || (to one-loop order, this has been noticed already in 
Refs.f^, 0). In the kinematical regime of interest, one finds || that the approximate 



fermion propagator and 3-point vertex are of the form: 

of i \ -z±{p) 7o =F P • 7 r i «± 7o ± P ■ 7 , OQ v 

S±(u;,p + q = — , T ± = . 28 

io - u±{p) - v±(p) • q 2 z± 2 

In these equations, the subscripts ± refer to the two positive-energy modes of the soft 
electronic, with energies w±(p) and residues z±(p), and v±(p) = dcu±(p)/dp. The main 
approximation used in the present paper remains meaningful for a soft quasiparticle whose 
momentum p ~ gT is much larger than the typical momenta q ^ g 2 T of the off-shell 
photons which are responsible for the leading IR singularity. By the same steps as before, 
we obtain the retarded propagator for the two fermionic modes ± in the form 

POO 

S±(w,p) = iz±(p) / dte^-" ±(p)+M?) A±(t), 
J o 

A ± (t) = Adutit), (29) 

with the function A(t) given by eq. (p3|). 

To conclude, we have shown that, in the weak coupling regime, quasiparticles ex- 
citations in the ultrarelativistic QED plasma are slowly damped. Their spectral den- 
sity retains the shape of a resonance strongly peaked around the perturbative mass-shell 
uo ~ v • p, with a typical width of order ~ g 2 T ln(l/p). No singularity is associated to this 
resonance: the retarded propagator is an entire function in the complex energy plane. 
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